High field magneto-transmission investigation of natural graphite 
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Magneto-transmission measurements in magnetic fields in the range B = 20 — 60 T have been 
performed to probe the H and Jf-point Landau level transitions in natural graphite. At the ff-point, 
two series of transitions, whose energy evolves as ^/B are observed. A reduced Slonczewski, Weiss 
and McClure (SWM) model with only two parameters to describe the intra-layer (70) and inter- 
layer (71) coupling correctly describes all observed transitions. Polarization resolved measurements 
confirm that the observed apparent splitting of the //-point transitions at high magnetic field cannot 
be attributed to an asymmetry of the Dirac cone. 

PACS numbers: 81.05.uf, 71.70.Di, 78.20.Ls, 78.30.-j 



o 



> 
o 



X 



Graphite consists of Bcrnal stacked sheets of hexag- 
onally arranged carbon atoms. The weak coupling be- 
tween the layers transforms the single graphene layer, 
which is a gapless semiconductor with a linear dispersion, 
into a semimetal with electron and hole puddles along 
the H — K — H edge of the hexagonal Brillouin zoneii 
In a magnetic field the electronic structure of graphite 
is accurately described by the Slonczewski, Weiss and 
McClure (SWM) band structure calculations?^ which 
require seven tight binding parameters 70, ..,75, A to 
define the interaction energy of the carbon atoms in 
the graphite lattice. The SWM model has been exten- 
sively verified using Shubnikov de Haas, de Haas van 
Alphen, thermopower and magneto-reflectance experi- 
mcnts4r— Carriers at the //-point behave as relativistic 
Dirac Fermions with a linear dispersion as in graphene. 
Magneto-absorption has been used to perform Landau 
level spectroscopy of carriers at the //-point {kz ~ 0.5) 
and /T-point (fc^ = 0) in both natural graphite and highly 
ordered pyrolytic graphite (HOPG)i^"— 

At the //-point, transitions with a characteristic \JnB 
magnetic field dependence of their energy are observed, 
which are identical to the transitions at the K and K' - 
points observed in graphene. For this reason, we refer 
to this series as "graphene-like" although we stress that 
here the series arises from the //-point transitions of per- 
fect bulk graphite. In addition, a second weaker series 
of transitions with a characteristic \JnB magnetic field 
dependence of their energy are observed. These transi- 
tions are absent in graphene, in fact they correspond to 
dipole forbidden transitions of the "graphene-like" series. 
However, this series correspond to dipole allowed transi- 
tions in graphite due to the complicated band structure 
at the //-pointi ^^d^d'^ We refer to these transitions as the 
"graphite-like" series, since they are absent in graphene. 
For the /T-point there is evidence of a splitting of the 
transitions which has been attributed to electron-hole 
asymmetryji^ 

Here we report magneto-optical absorption measure- 
ments to probe the evolution of the K and H point 
transitions in magnetic fields up to 60 T. This extends 
previous worki^rJ^ to higher magnetic fields and more 



importantly to higher energies. In particular, the use 
of near visible radiation facilitates the implementation of 
polarization resolved measurements. The observed trans- 
mission spectra are dominated by the Dirac-like series 
of transitions from the //-point. All the observed tran- 
sitions can be assigned, and the magnetic field evolu- 
tion reproduced, using a reduced SWM model with two 
tight binding parameters 70 and 71 . Polarization resolved 
measurements confirm that the observed splitting of the 
//-point "graphene-like" transitions is not linked to the 
asymmetry of the Dirac cone, which is anyway irrelevant 
at the //-point within the SWM model. Upon closer ex- 
amination, the splitting resembles rather an avoided level 
crossing, while nevertheless remaining unexplained. 

Thin samples for the transmission measurements were 
prepared by exfoliating natural graphite. The average 
thickness of the graphite layers remaining on the foil was 
estimated to be ~ 100 nm.— The measurements were 
performed up to 34 T at the dc resistive magnet labora- 
tory in Grenoble and up to 60 T at the pulsed magnetic 
field laboratory in Toulouse. For the absorption mea- 
surements a tungsten halogen lamp was used to provide 
broad spectrum in the visible and near infra-red range. 
The absorption was measured in the Faraday configura- 
tion in which fc, the wave propagation vector is parallel to 
the magnetic field, B. The c-axis of the graphite sample 
was parallel to magnetic field. A nitrogen cooled InGaAs 
photodiodc array coupled to a spectrometer collected the 
transmitted light from the sample in the spectral range 
850 — 1600 nm, i.e. energies of 0.8 — 1.5 eV. For the 
pulsed field measurements the exposure time was lim- 
ited to 2 ms in order to limit variations in the magnetic 
field during acquisition. Thirty spectra were taken dur- 
ing a 60 T shot so that in principle a complete magnetic 
field dependence can be acquired in a single shot. The 
magnetic field was systematically measured using a cal- 
ibrated pick-up coil. Since the absorption lines in this 
energy range are weak all the spectra were normalized 
by the zero field transmission to produce a differential 
transmission spectra. 

Typical differential magneto-absorption spectra mea- 
sured at r = 4.2 K for magnetic fields 48 - 58 T are 
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shown in Figure [Ua) . All spectra show a number of ab- 
sorption lines which can be assigned to dipole allowed 
transitions at the H and K points. The energetic po- 
sition of the observed absorption lines is plotted as a 
function of the magnetic field in Figure [TJb). In order to 
assign the transitions we first calculate the energy of the 
dipole allowed transitions (An = ±1) at the H and K- 
points using a greatly simplified SWM model with only 
two parameters 70 and A71 to describe the intra- and 
inter-layer coupling . ^^i^^i^^" — Here A = 2cos(7rfc2) and 
kz is the momentum perpendicular to the layers. This 
corresponds to treating graphite as a series of graphene 
bi- layers whose effective coupling depends on kz- The 
magneto-optical response is dominated by the singulari- 
ties in the joint density of initial and final states which 
occur at the iiT-point (A = 2) and iJ-point (A = 0). The 
energy spectrum of the Landau levels using the effective 
bilayer model is then given by, 

i?J±=±i=[(A7if + (2n + l)£2 

- V(A7i)^ + 2(2n + l)e2(A7i)2 + , (1) 



where e ~ c\/2ehB is the characteristic magnetic energy, 
c = '\/3eao7o/2?i is the Fermi velocity, ao = 0.246 nm 
is the lattice constant in the ab plane and ± labels the 
electron and hole Landau levels respectively. At the H- 
point, equation ^ reduces to the Landau level spectrum 
of graphene with E^^ = ±c\/2ehBn. 

The bi-layer model is expected to be almost exact at 
the iJ-point since the effect of trigonal warping (73) van- 
ishes and analytic expressions for the Landau levels can 
be easily obtained within the SWM model by diagonal- 
izing the Hamiltonian.— However, the situation is com- 
plicated by the presence of the Ei and E2 bands (see 
Fig 12), which are almost degenerate with E^ at the H- 
point (energy splitting A ~ —0.007 eV). In a magnetic 
field, neglecting the two exceptional E3 Landau levels 
(n = 0,-1), this gives rise to a second Landau level 
spectrum, E12 — E^±^ where n = 1,2,3,.. which is 
exactly degenerate with the i?3± ladder J^iiiii^ The Lan- 
dau level spectrum at the iJ-point is shown schematically 
in Fig. [5] where we indicate all possible dipole allowed 
^2(3) _^ ^3(2) transitions as an example. The graphene 

like transition E^_ — >■ E^\_ (labeled D) have the same 

energy as the £2^^"^ — eI'"^^ transitions which have a 
quantum number n which is lower by one. The circular 
polarization of the light required to excite each transition 
is indicated and we have adopted the convention that 
polarization corresponds to An = +1. The transition 
labeled d and d' are specific to graphite ("graphite-like" 
series). Transition d is the dipole allowed (|An| = 1) de- 
generate 'mixed' transitions E^_ — > E^ and i?| — >■ £^3+ 
which correspond to (are exactly degenerate with) dipole 
forbidden (An ~ 0) transitions of the graphene series. 
Transition d' shows dipole allowed (|An| ~ 1) degener- 
ate 'mixed' transitions E^_ — Ef and E2 E^_^_ which 
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FIG. 1: (Color online) (a) Typical differential magneto- 
transmission spectra of natural graphite measured at mag- 
netic fields in the range 48 - 58 T at T = 4.2 K. (b) Magnetic 
field dependence of the observed transitions assigned as fol- 
lows: //-point, graphene series (blue balls), graphite specific 
series (red triangles); /f-point (black stars). The lines are cal- 
culated energies of the dipole allowed //-point (solid, dotted, 
and, dot-dashed lines) and K point (dashed lines) transitions 
as described in the text. 



correspond to (are exactly degenerate with) dipole for- 
bidden transitions | An| = 2 of the graphene series. Note, 
that while we cannot exclude the presence in our sample 
of decoupled graphene layers, with transitions degenerate 
with the "graphene-like" series, the overwhelming contri- 
bution of graphite to the transmission is demonstrated by 
the observed strength of the "graphite-like" series. 

The energy of the dipole allowed optical transitions, 
calculated using Equation([T]) with the tight binding pa- 
rameters, 7o = 3.15 eV (c = 1.02 x 10^ m.s^^) and 
71 = 0.375 eV determined from magneto-absorption mea- 
surements at lower magnetic fieldsfi^ are plotted as a 
function of the magnetic field in Fig. [TJb) (solid and 
broken lines). The //-point transitions depend only on 
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FIG. 2: (Color online) (left) Band structure of graphite along 
the hexagonal H — K — H axis, (right) Dipole allowed transi- 
tions (An = 1) in a magnetic field at the iif-point correspond- 
ing the E2.- — ^ -E3+ (and E2 — >■ -Ei) "graphene-like" transi- 
tions (labeled D), and the 'mixed' E^- — >■ E2 and E2 — >■ -E3+ 
"graphite-like" transitions (labeled d and d'). The circular po- 
larization required to excite each transition is indicated using 
the convention that polarization corresponds to An = +1. 



the parameter 70 and evolve always as ^/B. The Ap- 
point transition depend also on the inter-layer coupling 
A71, and therefore, as can be seen from Equation([T]), 
evolve linearly at low energy (e <C A71) before increas- 
ing as ^fB at high energies (e ^ A71). At the i^-point 
A71 =0.75 eV so that we are in the intermediate regime 
where dependence is somewhere between linear and ^/B. 

The agreement between the reduced two-parameter 
SWM model and experiment in Figure [Ijb) is remark- 
able, especially taking into account that there are no ad- 
justable parameters. The iJ-point transitions iJ^'^' — > 
are labeled as in Figl^l Mainly if -point tran- 
sitions are observed, notably the "graphene-like" series 
i?3_ — ^ i?3+ and E2 — >■ El (thick blue solid lines la- 
beled with upper case letters) together with the weaker 
i?3_ — !• El and E2 — > -E3+ transitions (red dashed and 
dot-dashed lines labeled with lower case letters) . The K- 
point transitions, shown as black dotted lines are only 
observed directly at high magnetic fields. For com- 
pleteness, for the i^-point transitions, we have included 
phenomenologically the electron-hole asymmetry as sug- 
gested in Refs.dlill by using a different Fermi velocity 
Ce = 1.098 X lO'^ ms-i and cu = 0.942 x 10^ ms'^ for 
the electrons and holes respectively. These values arc 
slightly different from those used in Reffisl in order to 
have the same "average" value of 70 = 3.15 eV. While 



the electron-hole asymmetry was clearly seen in mea- 
surements at low magnetic fieldf^ the phenomenologi- 
cal asymmetry splitting introduced in Ref.[lB| decreases 
rapidly with increasing quantum number, and is probably 
too small to be seen in our high magnetic field data (the 
lowest energy A'-point transition seen is ri = 5, labeled 
i-5(-6) Lq(5) in Figure ma)). 

The "graphene-like" series unexpectedly shows what 
looks at first sight to be a splitting, which is puzzling 
since such a splitting is completely absent in magneto- 
transmission measurements on graphcnci^ This appar- 



ent splitting is clearly seen in the Agl"^^-* 
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transitions {n = 2,3,4) labeled D,E and F in FigH] 
However, a closer inspection of the magnetic field evolu- 
tion of the energy of the strong E^_ ' — !• E^]^ ' in Fig. [Ub) 
(transition D) indicates that the calculated transition fits 
better to the low energy feature at low fields {B < 50 T) 
before fitting better to the high energy feature at high 
fields {B > 54 T). This is suggestive of an avoided level 
crossing rather than a splitting. This hypothesis is sup- 
ported by the absorption spectra in Fig. [TJa), where it is 



clearly seen that the Agi^-* 



A3J|_ doublet (transition 



D) consists of a stronger low energy transition at low 
magnetic fields which switches to a stronger high energy 
feature at high fields, i.e. the two lines anti-cross. The 
origin of this behavior remains to be elucidated. How- 
ever, this cannot be due to inhomogeneity of the sample. 
A slightly different Fermi velocity for different regions 
would simply lead to an increased splitting with increas- 
ing magnetic field. 

Figure ^a,) shows differential absorption spectra mea- 
sured at i? = 58 T for different temperatures in the range 
4 — 300 K. A temperature of 100 K is already sufficient 
to suppress the apparent splitting of the E. 
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transition, although the transition itself, while weaken- 
ing slightly, remains clearly visible even at room tem- 
perature. The r = 4 K spectra have been fitted using 
a Lorentzian line shape of full width at half maximum 
(FWHM) of 28.5 meV for all transitions. The result of 
the fit (solid thin black line in Figl3Ja)) describes the 
data extremely well. The individual Lorentzians, for each 
transition are shown as dotted lines. Clearly, the broad- 
ening of the transitions is comparable to the energy sep- 
aration so that the absorption is only weakly modulated. 
Note that the disagreement between the data and the fit 
around 0.92 eV is probably a signature of the "missing" 
L_g(_7) — >L7(g) AT-point transition in Fig[l] which could 
not be assigned from the raw data. Keeping all other 
parameters constant, increasing the broadening of the 
Lorentzians produces a reasonable fit to the higher tem- 

2(3) 3f2) 

perature data, with the exception of the E^_ ' Ag^ ' 
transition. A reasonable fit to this transition at higher 
temperatures requires, in addition to a thermal broad- 
ening, that the amplitude of the two Lorentzian compo- 
nents be changed for which we see no physical justifi- 
cation. We therefore conclude that thermal broadening 
alone cannot explain the observed temperature depen- 
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FIG. 3: (Color online) (a) Differential magneto-transmission 
spectra of natural graphite measured at .B = 58 T for different 
temperatures. The thin black lines are a fit to the T — 4 K 
spectra assuming a Lorentzian line shape, (b) Polarization 
resolved spectra measured at .B = 59 T and T — 4.2 K for 
the transitions D, d and d' sketched schematically in Fig. [2] 
The polarization (cr^ or a^) has been arbitrarily assigned to 
a given magnetic field direction. 

dence of the E'g^^-' — > E^^"^ transition . 

^From a theoretical point of view a splitting of the 

-E'ai^'' — >■ E^+^ transition (or the degenerate E^'"'^^ — >■ 
2(1) 

E-^ ) transition) is not expected since even in the full 
SWM model. The effect of 73 (trigonal warping) vanishes 
at the 7?-point so that the energy levels are determined 
only by 79. The non vertical inter layer coupling term 
74, which induces electron- hole asymmetry, plays no role. 
This can be verified experimentally using the polarization 
resolved absorption measured at 59 T shown in Fig. ^ih) 
which focuses on transitions D, d and d! . There is no 
difference between (t+ and <j~ spectra for E^^^^ — >■ 
(transition D) confirming that the apparent doublet can- 
not under any circumstances be assigned to electron-hole 
asymmetry. In contrast, the 'mixed' i?3_ — Ei and 
E2 — > i?3+ transitions [d and d') show a marked de- 
pendence on the circular polarization with one of the 
transitions almost vanishing with either or a~ ex- 
citation. Using the polarization selection rules sketched 
in Fig 13] this can be explained provided one of the inter- 
band transitions (-E3- —J- Ei or E2 — £"3+) dominates. 



However, as in our experiment the sense of the circu- 
lar polarization has been arbitrarily assigned to a given 
magnetic field direction it is unfortunately not possible 
to know which transition prevails. 

While trigonal warping plays no role at the i/-point 
because 73 always enters the SWM Hamiltonian as 
73COs(7rfcz), close to the i7-point it can lead to mag- 
netic breakdown producing a splitting of levels in the 
Landau level structure which could possibly be observed 
in magneto-optical spectra at the iJ-point.— This orig- 
inates from an anti-crossing of Landau levels from the 
£"3 band with Landau levels from the Ei or E2 bands. 
The repulsion occurs due to the interaction caused by 73 
provided the Landau levels originate from the same sub- 
matrix (of the three possible) of the magnetic Hamilto- 
nian. In contrast to the iC-point, where trigonal warping 
induced magnetic breakdown occurs only at low mag- 
netic fields, close to the iJ-point magnetic breakdown 
takes place for all magnetic field strengths. An additional 
complication at very high magnetic fields [B « 70 T) is 
the predicted magnetic field induced transition of semi- 
metallic graphite to a zero gap semiconductor due to the 
crossing of the n = Landau level at the if-point and 
the n = —1 Landau level at the ff -point Further mea- 
surements at higher magnetic fields arc planned to clarify 
these issues. 

In conclusion, magneto-transmission measurements 
have been used to probe the H and if-point Landau 
level transitions in natural graphite. In the magnetic field 
range investigated, the spectra are dominated by transi- 
tions at the iJ-point. A "graphene-like" series together 
with a series of transitions exclusive to graphite are ob- 
served. We stress that both series arise from dipolc al- 
lowed transitions at the i7-point of perfect bulk graphite, 
and do not require the presence of decoupled graphene 
layers or decoupled bilayers in the sample. A reduced 
SWM model with only two parameter 70 and 71 correctly 
describes all observed transitions. Polarization resolved 
measurements (i) confirm that the apparent splitting of 
the "graphene-like" series at high magnetic field cannot 
be attributed to an asymmetry of the Dirac cone and (ii) 
suggest that the matrix elements connecting E^^ — > Ei 
and i?3_ — >■ i?2 are very different. 
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